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Summary. We study a discrete spatial model for invasive allele spread in which 
two alleles compete preemptively, initially only the "residents" (weaker competi- 
tors) being present. We find that the spread of the advantageous mutation is well 
described by homogeneous nucleation; in particular, in large systems the time- 
dependent global density of the resident allele is well approximated by Avrami's 
law. 



1.1 Introduction and Model 

The spatial and temporal characteristics of the spread of an advantageous 
mutation are fundamental questions in population dynamics. Fisher [1] and 
Kolmogorov et al. [2] first addressed these questions using the framework of 
a simple reaction-diffusion equation [3] . That work and many others focused 
on the velocity of the propagating front, which exists initially and separates 
the two spatial regions occupied separately by the two alleles. Both contin- 
uum and discrete spatial models have successfully tackled various aspects of 
these problems [3,4]. In this work we investigate how the advantageous allele 
emerges from "scratch" ; i.e., initially the region is fully dominated by the resi- 
dent allele and the advantageous allele is introduced by rare mutations. While 
the mutant allele has an individual-level advantage over the original one, the 
low probability of mutations, combined with a discrete spatial dynamics, can 
prevent the spread of the mutant for long times. Here we consider a model 
where the original "resident" and the competitively superior "invasive" allele 
compete for a common limiting resource preemptively [5-8]. 

The details of our model are as follows. We consider an LxL lattice with 
periodic boundary conditions. Each site can be empty or occupied by a single 
allele (either a resident or an invader). A lattice site represents the minimal 
level of locally available resource required to sustain an individual organism, 
hence the "excluded volume" constraint. We introduce the local occupation 
numbers at site x, rii(x) = 0, 1, i = 1, 2, representing the number of resident 
and invader alleles, respectively. By virtue of the excluded volume constraint, 
ni(x)n2(x) = 0. New individuals arise through local clonal propagation only. 
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That is, an individual occupying site x may reproduce if one or more neigh- 
boring sites are empty (here we consider nearest neighbor colonization only) . 
Competition for resources, hence space, is preemptive, therefore, an occupied 
site cannot be colonized by either allele until the current occupant's mor- 
tality leaves that site empty. Each individual may mutate and so carry the 
alternate allele; mutation is a two-way, recurrent process. 

We performed dynamic Monte Carlo simulations to study the above 
model. Our time unit is one Monte Carlo step per site (MCSS) during which 
L 2 sites are chosen randomly. If a site is empty, it may be colonized by indi- 
viduals of allele i occupying neighboring sites, at the rate o^r^x); at is the 
individual-level colonization rate and ryi(x) = (1/4) X)x'enn(x) n i( x ') i s the 
density of allele i around site x [nn(x) is the set of nearest neighbors of site 
x] . If a site is occupied by an individual, it can die at rate p, (regardless of the 
allele) or mutate to the other allele at rate 4>. We can summarize the local 
transition rules for an arbitrary site x as 

O^'l, 0" 2 -H x) 2, 1-^0, 2-^0, 1^2, (1.1) 

where 0,1,2 indicates whether the site is empty, occupied by an individual 
with the resident, or an individual with the invader allele, respectively. 

In the simulations, we initialized the system fully occupied by the resident 
allele (ni(x)=l for all x). We are interested in the parameter region where 
4> <C p < a\ < a.2, so that mutation is a rare process, but the invader 
allele has a reproductive-effort advantage. Then, due to mortality, the system 
quickly relaxes (much too fast for mutation to play a role) to the "quasi- 
cquilibrium" state where the resident's population is balanced by its own 
clonal propagation and mortality rates (in the near absence of invaders). 
Throughout the simulations, we track the time-dependent global densities 
of the two alleles, pi(t) = (1/L 2 ) J2 X n i( x , t). We define the lifetime r of 
the resident allele as the first passage time of pi (t) to one-half of its quasi- 
cquilibrium value p\. 

1.2 Single-Cluster and Multi-Cluster Spread 

As a result of the rare mutations, individuals with the invasive allele occasion- 
ally appear in the population. An invader lacking access to nearby resources 
may die without propagating. If a site opens in the local neighborhood (re- 
source becomes available), the invader may colonize it. However, the empty 
site is likely surrounded by more than one resident. The resident's greater lo- 
cal density can compensate for its lower individual-level colonization rate, so 
the resident has the better chance of colonizing an empty site. Consequently, 
one expects small clusters of the invading allele to shrink and disappear. Resi- 
dents, although weaker competitors, can prevail for some time, since preemp- 
tive competition imposes a strong constraint on the growth of the invaders. 
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Fig. 1.1. (a) Cumulative probability distributions P no t{t) for L—32, Qi=0.50, 
Q2=0.70, and fi— 0.20 for three different values of the mutation rate <f> (in increas- 
ing order from the top), (b) Average nucleation time (in units of MCSS) vs. the 
mutation rate for two different values of n [ai, ai, and L are the same as for (a)]. 
The inset shows the same on log-log scales. The straight solid line corresponds to 
a slope —1, indicating (t n ) ~ in the single-cluster regime. 

Individuals with the advantageous allele can succeed only if they generate a 
cluster large enough that it statistically tends to grow at its periphery. 

Snapshots of configurations and preliminary studies [9] confirm the exis- 
tence of a critical cluster size, beyond which the spread of the invading allele 
becomes statistically favorable. Further, they also show strongly clustered 
growth of the invading allele. For a given set of parameters, there exists a 
length scale R Q , the typical spatial separation of invading clusters; for L<S^R 
the invasion almost always occurs through the spread of a single invading clus- 
ter [single-cluster (SC) invasion], while for L^>R Q the invasion is the result of 
many invading clusters [multi-cluster (MC) invasion]. Conversely, fixing the 
linear system size L and other parameters (except the mutation rate) there 
is a characteristic value of (f> (now controlling R a ), such that for sufficiently 
low values of </>, MC invasion by the advantageous allele crosses over to the 
SC pattern. 

The above picture suggests that we can apply the framework of homo- 
geneous nucleation and growth [10-12] to describe the spatial and temporal 
characteristics of the spread of the invasive allele. This framework has suc- 
cessfully described analogous dynamic phenomena in ferromagnetic [13-15] 
and ferroelectric materials [16,17], flame propagation in slow combustion [18], 
chemical reactions [19], and other ecological systems [9,20,21]. While local 
mutation is a Poisson process, lacking a Hamiltonian or an effective free en- 
ergy for the model, it is not known a priori whether the nucleation of a "su- 
percritical" cluster will also be Poisson. To this end, in the SC regime, we con- 
structed cumulative probability distributions for the lifetime of the resident 
allele P no t(t), i.e., the probability that the global density of the resident has 
not crossed below p\/2 by time t. We found that these distributions are indeed 
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Fig. 1.2. (a) Time-series of the global density of the resident allele pi(t) for 
L—l, 000, ai=0.50, a2=0.70, and /i=0.20, for three different values of the mutation 
rate (j> (in increasing order from the top). The solid curves represent Avramis's law, 
Eq. (1.2). The inset shows the pi(t) vs. t 3 on log-linear scales, (b) Average lifetime 
(in units of MCSS) vs. the mutation rate on log-log scales for the same parameter 
values. The straight dashed line is the best fit power-law indicating (r) ~ (f>~°' S04 . 



exponentials (indicating that the nucleation of a successful invading cluster 
is a Poisson process): P no t(t) = 1 for t<t g and P no t{t) = cxp[— (t — t g )/(t n )] 
for t>t g . Here (t n ) is the average nucleation time and t g is the close-to- 
deterministic growth time until the advantageous mutation dominates half 
the system. We show results for a fixed (sufficiently small) system size for 
three mutation rates in Fig. 1.1(a). From the slopes of the exponentials we 
obtained the average nucleation times [Fig. 1.1(b)], hence the ^-dependence 
of the nucleation rate per unit volume /(</>)• Since (t n ) — [L 2 /(0)]~ 1 , we 
have I(4>) ~ (tn)^ 1 ~ </>• In the SC regime, the invasive spread is inher- 
ently stochastic; it is initiated and completed by the first randomly nucle- 
ated successful cluster of the advantageous allele. For very low values of <j>, 
the lifetime is dominated by the very large average nucleation times, hence 
(r) = (t n )+t g « (U^" 1 . 

In the MC regime the invasion processes becomes self-averaging and the 
global densities approach deterministic functions in the limit of L^oo. At the 
same time, (r) approaches a system-size independent limit. For large systems 
we applied the KJMA theory [10-12] (or Avrami's law) to predict the density 
of the resident alelle, 

Pl {t) ~p* e -M2)(V(-» 3 . (i.2) 

Our results in Fig 1.2(a) show that it is, indeed, a very good approximation. 
Assuming that the spreading velocity of the invading clusters is constant, 
KJMA theory predicts that (r) ~ [/((/>)]~ 1/3 - </>~ 1/3 in the MC regime. The 
measured exponent, —0.304, is not too far off [Fig. 1.2(b)], but it indicates 
that the assumption of a constant spreading velocity (possibly as a result of 
the nontrivial surface properties of the clusters) may break down. 
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I. 3 Summary and Outlook 

We studied the spread of an advantageous mutant in a two- allele population 
where rare mutations introduce the favored allele. We found that nuclcation 
theory, in particular Avrami's law, describes this phenomenon very well. Sys- 
tematic studies of the critical cluster size and the cluster-size dependence of 
the spreading velocity are under way. The structure of the spreading clusters, 
in particular, the roughness of their surface, is expected to play an important 
role in the latter. 

References 

1. R.A. Fisher: Annals of Eugenics 7, 355 (1937) 

2. A.N. Kolmogorov, I. Petrovsky, N. Piskounov: Moscow Univ. Bull. Math. 1, 1 
(1937) 

3. J.D. Murray: Mathematical Biology I and II, 3rd edition (Springer-Verlag, New 
York, 2002, 2003) 

4. D. ben-Avraham, S. Havlin: Diffusion and Reactions in Fractals and Disor- 
dered Systems (Cambridge University Press, Cambridge, UK, 2000) 

5. J.B. Shurin, P. Amarasekare, J.M. Chase, R.D. Holt, M.F. Hoopes, M.A. Lei- 
bold: J. Theor. Biol. 227, 359 (2004) 

6. P. Amarasekare: Ecol. Lett. 6, 1109 (2003) 

7. D.W. Yu, H.B. Wilson: Am. Nat. 158, 49 (2001) 

8. D.E. Taneyhill: Ecol. Monogr. 70, 495 (2000) 

9. L. O'Malley, A. Allstadt, G. Korniss, T. Caraco: in Fluctuations and Noise in 
Biological, Biophysical, and Biomedical Systems III, edited by N.G. Stocks, D. 
Abbott, and R.P. Morse, Proceedings of SPIE Vol. 5841 (SPIE, Bellingham, 
WA, 2005, in press) 

10. A.N. Kolmogorov: Bull. Acad. Sci. USSR, Phys. Ser. 1, 355 (1937) 

II. W. A. Johnson and R. F. Mehl: Trans. Am. Inst. Mining and Metallurgical 
Engineers 135, 416 (1939) 

12. M. Avrami: J. Chem. Phys. 7, 1103 (1939); J. Chem. Phys. 8, 212 (1940); J. 
Chem. Phys. 9, 177 (1941) 

13. P.A. Rikvold, H. Tomita, S. Miyashita, S.W. Sides: Phys. Rev. E 49, 5080 
(1994) 

14. H.L. Richards, S.W. Sides, M.A. Novotny, P.A. Rikvold: J. Magnet. Mater. 
150, 37 (1995) 

15. R.A. Ramos, P.A. Rikvold, M.A. Novotny: Phys. Rev. B 59, 9053 (1999) 

16. Y. Ishibashi, Y. Takagi: J. Phys. Soc. Jpn. 31, 506 (1971) 

17. H.M. Duiker, P.D. Beale: Phys. Rev. B 41, 490 (1990) 

18. M. Karttunen, N. Provatas, T. Ala-Nissila, M. Grant: J. Stat. Phys. 90, 1401 
(1998) 

19. E. Machado, G. M. Buendia, P. A. Rikvold: Phys. Rev. E (in press, 2005); 
arXiv: cond-mat/0411335 

20. G. Korniss, T. Caraco: J. Theor. Biol. 233, 137 (2005) 

21. A. Gandhi, S. Levin, S. Orszag: J. Theor. Biol. 200, 121 (1999) 

22. This research is supported in part by the US NSF through Grant Nos. 
DMR-0113049, DMR-0426488, DEB-0342689 and by the Research Corpora- 
tion through Grant No. RI0761. 



